The Polarizable Continuum Model Goes Viral! Extensible, Modular and Sustainable Development of Quantum Mechanical Continuum Solvation Models

Synergistic theoretical and experimental approaches to challenging chemical problems have become more and more widespread, due to the availability of efficient and accurate ab initio quantum chemical models. Limitations to such an approach do, however, still exist. The vast majority of chemical phenomena happens in complex environments, where the molecule of interest can interact with a large number of other moieties, solvent molecules or residues in a protein. These systems represent an ongoing challenge to our modelling capabilities, especially when high accuracy is required for the prediction of exotic and novel molecular properties. How to achieve the insight needed to understand and predict the physics and chemistry of such complex systems is still an open question. I will present our efforts in answering this question based on the development of the polarizable continuum model for solvation. While the solute is described by a quantum mechanical method, the surrounding environment is replaced by a structureless continuum dielectric. The mutual polarization of the solute-environment system is described by classical electrostatics. Despite its inherent simplifications, the model contains the basic mathematical features of more refined explicit quantum/classical polarizable models. Leveraging this fundamental similarity, we show how the inclusion of environment effects for relativistic and nonrelativistic quantum mechanical Hamiltonians, arbitrary order response properties and high-level electron correlation methods can be transparently derived and implemented. The computer implementation of the polarizable continuum model is central to the work presented in this dissertation. The quantum chemistry software ecosystem suffers from a growing complexity. Modular programming offers an extensible, flexible and sustainable paradigm to implement new features with reduced effort. PCMSolver, our open-source application programming interface, can provide continuum solvation functionality to any quantum chemistry software: continuum solvation goes viral. Our strategy affords simpler programming workflows, more thorough testing and lower overall code complexity. As examples of the flexibility of our implementation approach, we present results for the continuum modelling of non homogeneous environments and how wavelet-based numerical methods greatly outperform the accuracy of traditional methods usually employed in continuum solvation models

Diagrammatic Coupled Cluster Monte Carlo

We propose a modified coupled cluster Monte Carlo algorithm that stochastically samples connected terms within the truncated Baker--Campbell--Hausdorff expansion of the similarity transformed Hamiltonian by construction of coupled cluster diagrams on the fly. Our new approach -- diagCCMC -- allows propagation to be performed using only the connected components of the similarity-transformed Hamiltonian, greatly reducing the memory cost associated with the stochastic solution of the coupled cluster equations. We show that for perfectly local, noninteracting systems, diagCCMC is able to represent the coupled cluster wavefunction with a memory cost that scales linearly with system size. The favorable memory cost is observed with the only assumption of fixed stochastic granularity and is valid for arbitrary levels of coupled cluster theory. Significant reduction in memory cost is also shown to smoothly appear with dissociation of a finite chain of helium atoms. This approach is also shown not to break down in the presence of strong correlation through the example of a stretched nitrogen molecule. Our novel methodology moves the theoretical basis of coupled cluster Monte Carlo closer to deterministic approaches.Comment: 31 pages, 6 figure

Fully Polarizable QM/Fluctuating Charge Approach to Two-Photon Absorption of Aqueous Solutions

We present the extension of the quantum/classical polarizable fluctuating charge model to the calculation of single residues of quadratic response functions, as required for the computational modeling of two-photon absorption cross-sections. By virtue of a variational formulation of the quantum/classical polarizable coupling, we are able to exploit an atomic orbital-based quasienergy formalism to derive the additional coupling terms in the response equations. Our formalism can be extended to the calculation of arbitrary order response functions and their residues. The approach has been applied to the challenging problem of one- and two-photon spectra of rhodamine 6G (R6G) in aqueous solution. Solvent effects on one- and two-photon spectra of R6G in aqueous solution have been analyzed by considering three different approaches, from a continuum (QM/PCM) to two QM/MM models (non-polarizable QM/TIP3P and polarizable QM/FQ). Both QM/TIP3P and QM/FQ simulated OPA and TPA spectra show that the inclusion of discrete water solvent molecules is essential to increase the agreement between theory and experiment. QM/FQ has been shown to give the best agreement with experiments

Diagrammatic Coupled Cluster Monte Carlo.

We propose a modified coupled cluster Monte Carlo algorithm that stochastically samples connected terms within the truncated Baker-Campbell-Hausdorff expansion of the similarity-transformed Hamiltonian by construction of coupled cluster diagrams on the fly. Our new approach-diagCCMC-allows propagation to be performed using only the connected components of the similarity-transformed Hamiltonian, greatly reducing the memory cost associated with the stochastic solution of the coupled cluster equations. We show that for perfectly local, noninteracting systems diagCCMC is able to represent the coupled cluster wavefunction with a memory cost that scales linearly with system size. The favorable memory cost is observed with the only assumption of fixed stochastic granularity and is valid for arbitrary levels of coupled cluster theory. Significant reduction in memory cost is also shown to smoothly appear with dissociation of a finite chain of helium atoms. This approach is also shown not to break down in the presence of strong correlation through the example of a stretched nitrogen molecule. Our novel methodology moves the theoretical basis of coupled cluster Monte Carlo closer to deterministic approaches.Sims Fun

The Distinct Cases of Kosovo and South Ossetia: Deciding the Question of Independence on the Merits and International Law

The cases of Kosovo and South Ossetia provide two opportunities for the community of nations to reconcile the rights of Serbia and Georgia with the rights of the peoples within their borders. Instead of doing so, other countries used these cases for ideological and political posturing and continued to do so as of the end of 2008. This Note will review Kosovo and South Ossetia and attempt to take the politics out of an inherently political decision-whether or not to recognize them as independent states. Part I of this Note will review how the United Nations has approached the issues of autonomy, borders, secession, and minorities. It will also provide background on the doctrines of self-determination and the rights of states. Part II will set out the facts-as they are generally accepted-of Kosovo and South Ossetia and consider the differing opinions regarding independence of the two areas. Part III will review those opinions in light of the goals and doctrines of the United Nations and how it has approached these issues in the past. The Note concludes that recognition of independence of the states should be driven by legal arguments rather than political or social alliances and that the goals of the United Nations are best served by being steadfast in promoting multi-ethnicity within political entities; Kosovo and South Ossetia should be recognized as autonomous regions within Serbia and Georgia, respectively. The Note also concludes that domestic as well as international legal structures must be in place to make autonomy workable

Cavity-free continuum solvation: implementation and parametrization in a multiwavelet framework

We present a multiwavelet-based implementation of a quantum/classical polarizable continuum model. The solvent model uses a diffuse solute-solvent boundary and a position-dependent permittivity, lifting the sharp-boundary assumption underlying many existing continuum solvation models. We are able to include both surface and volume polarization effects in the quantum/classical coupling, with guaranteed precision, due to the adaptive refinement strategies of our multiwavelet implementation. The model can account for complex solvent environments and does not need a posteriori corrections for volume polarization effects. We validate our results against a sharp-boundary continuum model and find very good correlation of the polarization energies computed for the Minnesota solvation database

The HANDE-QMC Project: Open-Source Stochastic Quantum Chemistry from the Ground State Up.

Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included

The DIRAC code for relativistic molecular calculations

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree?Fock, Kohn?Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.Fil: Saue, Trond. Université Paul Sabatier; Francia. Centre National de la Recherche Scientifique; FranciaFil: Bast, Radovan. Uit The Arctic University Of Norway; NoruegaFil: Gomes, André Severo Pereira. University Of Lille.; Francia. Centre National de la Recherche Scientifique; FranciaFil: Jensen, Hans Jorgen Aa.. University of Southern Denmark; DinamarcaFil: Visscher, Lucas. Vrije Universiteit Amsterdam; Países BajosFil: Aucar, Ignacio Agustín. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Modelado e Innovación Tecnológica. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas Naturales y Agrimensura. Instituto de Modelado e Innovación Tecnológica; Argentina. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas y Naturales y Agrimensura. Departamento de Física; ArgentinaFil: Di Remigio, Roberto. Uit The Arctic University of Norway; NoruegaFil: Dyall, Kenneth G.. Dirac Solutions; Estados UnidosFil: Eliav, Ephraim. Universitat Tel Aviv.; IsraelFil: Fasshauer, Elke. Aarhus University. Department of Bioscience; DinamarcaFil: Fleig, Timo. Université Paul Sabatier; Francia. Centre National de la Recherche Scientifique; FranciaFil: Halbert, Loïc. Centre National de la Recherche Scientifique; Francia. University Of Lille.; FranciaFil: Hedegård, Erik Donovan. Lund University; SueciaFil: Helmich-Paris, Benjamin. Max-planck-institut Für Kohlenforschung; AlemaniaFil: Ilias, Miroslav. Matej Bel University; EslovaquiaFil: Jacob, Christoph R.. Technische Universität Braunschweig; AlemaniaFil: Knecht, Stefan. Eth Zürich, Laboratorium Für Physikalische Chemie; SuizaFil: Laerdahl, Jon K.. Oslo University Hospital; NoruegaFil: Vidal, Marta L.. Department Of Chemistry; DinamarcaFil: Nayak, Malaya K.. Bhabha Atomic Research Centre; IndiaFil: Olejniczak, Malgorzata. University Of Warsaw; PoloniaFil: Olsen, Jógvan Magnus Haugaard. Uit The Arctic University Of Norway; NoruegaFil: Pernpointner, Markus. Kybeidos Gmbh; AlemaniaFil: Senjean, Bruno. Universiteit Leiden; Países BajosFil: Shee, Avijit. Department Of Chemistry; Estados UnidosFil: Sunaga, Ayaki. Tokyo Metropolitan University; JapónFil: van Stralen, Joost N. P.. Vrije Universiteit Amsterdam; Países Bajo

Proceedings of the SC15 Workshop on Producing High Performance and Sustainable Software for Molecular Simulation

These are the proceedings of the SC15 Workshop on Producing High Performance and Sustainable Software for Molecular Simulation, which was organised as part of the APES project. Please see the Introduction and Table of Contents for more details